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ABSTRACT 

Dark matter haloes in ACDM simulations grow by mergers with other haloes as well as 
accretion of "diffuse" non-halo material. We quantify the mass growth rates via these 
two processes, Mmor and Mdif, and their respective dependence on the local halo en- 
vironment using the ~ 500,000 haloes of mass ~ 10^^ to IO^^Mq in the Millennium 
simulation. Adopting a local mass density parameter as a measure of halo environ- 
ment, we find the two rates to show strong but opposite environmental dependence, 
with mergers playing an increasingly important role for halo growths in overdense 
regions while diffuse accretion dominating the growth in the voids. For galaxy-scale 
haloes, these two opposite correlations largely cancel out, but a weak environmental 
dependence remains that results in a slightly lower mean total growth rate, and hence 
an earlier mean formation redshift, for haloes in denser environments. The mean for- 
mation redshift of cluster-mass haloes, on the other hand, shows no correlation with 
halo environment. The origin of the positive correlation of Mmer with local density 
can be traced to the sourrounding mass reservoir outside the virial radii of the haloes, 
where more than 80% of the mass is in the form of resolved haloes for haloes residing 
in densest regions, while this fraction drops to ~ 20% in the voids. The negative cor- 
relation of Mdif with local density, however, is not explained by the available diffuse 
mass in the reservoir outside of haloes, which is in fact larger in denser regions. The 
non-halo component may therefore be partially comprised of truly diffuse dark mat- 
ter particles that are dynamically hotter due to tidal stripping and are accreted at a 
suppressed rate in denser regions. We also discuss the implications of these results for 
how to modify the analytic Extended Press-Schechter model of halo growths, which 
in its original form does not predict environmental dependence. 



1 INTRODUCTION 



The bottom-up growth of structure is a liallmark of hierar- 
chical cosmological models such as the A cold dark matter 
(ACDM) model. In these universes, dark matter haloes of 
lower mass are expected to form earlier, on average, than 
massive haloes. In observations, theoretical studies, and 
semi-analytical modelling of galaxy formation (see |Baugh| 
12006 for a review), the mass of a halo is the key variable 



of these studies, halo environment is characterised via the 
halo bias parameter, which is determined from the relative 
clustering strengths of haloes to the underlying dark matter 
distribution. 



upon which many properties of galaxies and their host haloes 
depend, e.g., formation redshift, galaxy occupation number, 
colour, morphology, star formation rate, and stellar feedback 
processes. 

Recent work has shown that in addition to the mass, 
various properties of halo formation and evolution also de- 
pend on the environment within which the haloes reside. For 
instance, at a fixed halo mass, older haloes have been found 
to be more clustered than younger haloes, and the correla- 
tion between clustering strength and formation redshift is 



In Part I of this series ( |Fakhouri fc Ma|2009| henceforth 
FM09), we have instead chosen to use local overdensities as 
a more direct measure of halo environment. The focus of 
the study was to quantify the environmental dependence of 
the merger rate of haloes. We compared a number of local 
overdensity variables (both including and excluding the halo 
mass itself), some of which had been used in earlier stud- 
ies ([Lemson fc Kauffmann|[T999j [Harker et al.||2006| |Wang| 



et al. 2007 Maulbetsch et al. 2007, Hahn et al. 20081. The 
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key finding in FM09 was that halo-halo mergers occur more 
frequently in denser regions than in voids, and that this en- 
vironmental dependence is similar regardless of the merger 
mass ratio (e.g., minor vs major) or the descendant halo 
mass (galaxy- vs cluster-sized): we found mergers to occur 
about 2.5 times more frequently in the densest regions than 
in the emptiest regions. We provided an analytical formula 
as a function of local density to approximate this environ- 
mental trend. This expression can be used with the fit for 
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the global mean rate of Fakhouri & Ma ( 2008 1 (henceforth 



tion of the surrounding mass reservoir outside of the virial 



FM08) to predict halo merger rates as a function of descen- 
dant mass, progenitor mass ratio, redshift, and environment 
over a wide range of parameter space. 

In this paper, we build on our earlier merger rate studies 
by investigating the mass growth of haloes, and its environ- 
mental dependence, via two sources: growth from mergers 
with other haloes (a quantity closely related to the results 
of FM09), and growth from accretion of non-halo material, 
which we will refer to as "diffuse" accretion. Due to the finite 
resolution of the simulations, we expect a portion of this dif- 
fuse component to be comprised of unresolved haloes, which, 
in a higher-resolution simulation, should largely follow the 
merger physics and scaling laws of their higher mass counter- 
parts in our earlier studies. The diffuse non-halo component, 
however, can in principal also contain truly diffuse dark mat- 
ter particles that either were tidally stripped from existing 
haloes or were never gravitationally bound to any haloes. 

As reported below, we find that diffuse accretion plays 
an important role in contributing to halo mass growth in 
the Millennium simulation. Moreover, we will show that this 
growth component correlates with the local halo environ- 
ment in an opposite way from the component due to merg- 
ers, with diffuse accretion playing an increasingly important 
role in halo growth in the voids, and mergers playing a more 
important role in the densest regions. This difference sug- 
gests that the diffuse component is not simply an extension 
of the resolved haloes down to lower masses, but rather that 
there is an intrinsic difference between the two components 
that results in the opposite environmental trends. An impli- 
cation of this result is that Milky- Way size galaxies in voids 
and those that reside near massive clusters may have statisti- 
cally distinguishable formation history, where the galaxies in 
voids acquire their baryons more quiescently via diffuse ac- 
cretion, while those in dense regions assemble their baryons 
mainly via mergers. Such environmental effect can show up 
in galaxy properties such as the star formation rates, colors, 
and morphologies. 

The environmental dependence of halo growths re- 
ported in this paper also has far-reaching implications for 
the much-used analytic theories for halo growth such as the 
Extended Press-Schechter (EPS) and excursion set models 
UPress fc Schechter|[T974} [Bond et al.|p9T] [Lacey fc Cole| 
19931. These models assume that all dark matter particles 



radii of the haloes (§ 4.4). \ 4.5 discusses a test that we have 



reside in haloes, and halo growths depend only on mass and 
not environment. As we will elaborate on below, both as- 
sumptions are too simplistic and must be modified to ac- 
count for the results from numerical simulations. 

This paper is organised as follows. §2 . 1 1 summarises the 
various definitions of halo mass and environment used in 
our analysis, as well as the means by which we extract halo 
merger trees from the public data in the Millennium simu- 
lation. In Sec. 3 we discuss how the two mass growth rates 
due to mergers and diffuse accretion, Mmor and A/dif, are 
defined and computed. The distributions of the rates for the 
~ 500, 000 haloes at 2 = and their redshift evolution are 
presented. The environmental dependence of halo growths 
is analysed in Sec. 4. The correlations of four halo proper- 
ties with the local density parameter are investigated: the 
mass growth rates Mmer and Mdif (§ 4.1), the fraction of 
a halo's final mass gained via mergers vs. diffuse accretion 
(§ 4.2), the formation redshift 2/ (§ 4.3), and the composi- 



performed to verify that the majority of haloes reside in a 
similar environmental region (e.g. overdense or underdense) 
throughout their lifetimes. In Sec. 5 we investigate further 
the nature of the "diffuse" component by varying the mass 
threshold used to define Mmor vs Afdif- We then discuss the 
implications of our results for the analytic EPS model of 
halo growth, which is entirely independent of environment 
in its basic form and therefore must be modified to account 
for the various environmental trends reported in § 4. 



2 MERGER TREES AND HALO 

ENVIRONMENT IN THE MILLENNIUM 
SIMULATION 

2.1 Merger Trees of the Millennium Haloes 



The Millennium simulation ( Springel et al. 2005 I assumes 
a ACDM model with f7„ = 0.25, VLi = 0.045, = 0.73, 
h = 0.73 and a spectral index of n = 1 for the primordial 
density perturbations with normalisation erg = 0.9 ( Springel 



et al. 20051. The dark matter N-body simulation followed 



the trajectories of 2160"^ particles of mass 1.2 x m a 

(685 Mpc)^ box from redshift z = 127 to z = 0. A friends-of- 
friends group finder with a linking length of fe = 0.2 is used 
to identify ~ 2 x 10^ dark matter haloes in the simulation 
down to a mass resolution of 40 particles (~ 4.7 x IO^^Mq). 
Each FOF halo thus identified is further broken into con- 
stituent subhaloes, each with at least 20 particles, by the 
SUBFIND algorithm that identifies gravitationally bound 



substructures within the host FOF halo (Springel et al. 
|2001) . 

Even though the Millennium public database provides 
a catalogue of FOF haloes at each output, it does not give 
the merger trees for these haloes. Instead, it provides merger 
trees for the subhaloes, which are constructed by connecting 
the subhaloes across the 64 available redshift outputs. Dur- 
ing this construction, a decision must be made about the 
ancestral relations of the subhaloes since the particles in a 
given subhalo may go into more than a single subhalo in the 
subsequent output. In this case, a subhalo is chosen to be 
the descendent of a progenitor subhalo at an earlier output 
if it hosts the largest number of bound particles in the pro- 
genitor subhalo. The resulting merger tree of subhaloes can 
then be processed further to construct the merger tree of 
the FOF haloes. This construction is non-trivial due to the 
fragmentation of FOF haloes; this is discussed at length in 
FM08 and FM09. 

In FM08 and FM09, we proposed a variety of post- 
processing algorithms to handle FOF fragmentation. Three 
methods were compared in FM08: snip, stitch-3, and stitch- 
00. The snipping method severs the ancestral relationship 
between halo fragments and their progenitors. This is the 
method commonly used in the literature and suffers from in- 
flated merger rates due to the aberrant remerger of snipped 
fragments. The two stitching algorithms prevent halo frag- 
mentation by "stitching" halo fragments together such that 
each FOF halo in the simulation has exactly one descendant. 
Stitch-oo performs this procedure whenever fragmentation 
occurs, whereas stitch-3 only reincorporates halo fragments 



© 0000 RAS, MNRAS 000, [T]-[T4l 



Environmental Dependence of Halo Growth 3 



Mass Percentile (Masses in units of IQ-'^^AfQ) 



Redshift 


0-40% 


40-70% 


70-90% 


90-99% 


99-100% 


2 = (Az = 0.06, At = 0.83 Gyr) 


1.2 - 2.1 


2.1 - 4.5 


4.5 - 14 


14- 110 


> 110 


z = 0.51 (Az = 0.06, At = 0.38 Gyr) 


1.2 - 2.0 


2.0 - 4.1 


4.1 - 12 


12 - 74 


> 74 


z = 1.08 (Az = 0.09, At = 0.34 Gyr) 


1.2 - 1.9 


1.9 - 3.7 


3.7-9.5 


9.5 - 48 


> 48 


z = 2.07 (Az = 0.17, At = 0.25 Gyr) 


1.2 - 1.8 


1.8 - 3.0 


3.0 - 6.5 


6.5 - 24 


> 24 



Table 1. Mass bins used in this paper for redshifts 0,0.51,1.08, and 2.07. The bins are computed assuming fixed mass-perccntile bins 
(header row) and the mass boundaries are computed using the prescribed percentile. The high mass bins extend out to 5.2 X 10^^,3 X 
10^^,1.3 X 10^^, and 4.4 X IO^^Mq for z = 0,0.51,1.08,2.07 respectively, though these exceptionally high mass objects are outliers and 
not particularly representative of the high mass range. 



that are destined to remerge within 3 outputs of the frag- 
mentation event. Both stitching algorithms lower the minor 
merger rate as they prevent spurious remergers of halo frag- 
ments. Naturally, this reduction is strongest for stitch-oo. 

In FM09 we showed that the rate of halo fragmentation 
was a strong function of halo environment, with haloes re- 
siding in overdense regions undergoing fragmentations over 
three times as often as haloes in underdense regions. We 
showed that the choice of post-processing algorithm was 
important in determining the environmental dependence of 
halo merger rates in dense environments and introduced two 
new algorithms: split-3 and split-oo. These algorithms re- 
move halo fragmentations by splitting progenitor haloes that 
fragment into multiple haloes to ensure that each halo has 
exactly one descendant. As with stitch-3, split-3 only splits 
progenitors that are split within the past 3 outputs (looking 
backwards towards high z) of the fragmentation event. Split- 
oo always splits fragmenting haloes. See also |Genel et al.| 
1 2008 1 who introduced a similar algorithm. 



We showed in FM09 that split-oo suffered from an "un- 
zipping" effect as halo fragmentations propagated up the 
tree, breaking the self-similarity in the merger rate, B/n, 
observed in FM08. Split-3, however, yielded very similar 
merger rates (within ~ 10%) as stitch-3 (see Fig. 10 of 
FM09). 

We have tested the various quantities presented in this 
paper using all five post-processing algorithms and have in- 
deed found the differences not to be significant enough to 
warrant presentating plots for all methods. We will, how- 
ever, briefly discuss any relevant differences in the text. The 
results presented throughout this paper will use the stitch-3 
algorithm used in FM08 and FM09. 



2.2 Halo Environment 

To measure each halo's local environment, we will primarily 
use the overdensity variable ^t-pop presented in FM09. This 
quantity computes the overdensity within a sphere of radius 
7 Mpc, centred on each halo in the simulation, by adding 
up the mass contributions made by the particles within the 
sphere and subtracting the FOF mass of the central halo. 
The database provides the dark matter density on a 256'' 
grid but not the full particle distribution itself. We have 
mapped this grid, which is given along a Peano-Hilbert space 
filling curve, to spatial coordinates and interpolated it to 
obtain 5 at the position of each halo's centre. 

We will also compare some of the results in this pa- 



per to the alternative environmental variable Sr discussed 
in FM09, where 5? is the overdensity surrounding a halo 
within a sphere of radius 7 Mpc (including the halo mass 
itself). Fig. 1 of FM09 showed that St was strongly corre- 
lated with mass for haloes above 10^''A/q when the central 
object's contribution to S began to dominate at these mass 
scales. This correlation was absent when the variable ^t-pop 
was used. We will therefore use (Jy-pop as the environmental 
measure for most of the analyses below. 

2.3 Halo Mass Definitions 

The Millennium database provides a number of mass mea- 
surements for each identified FOF halo. Two simple mass 
measures are Mpop, the total mass of the particles con- 
nected to an FOF by the group finder, and Msh, the sum 
of the masses of the subhaloes that constitute the FOF halo. 
The latter mass includes only particles that are bound by 
the SUBFIND algorithm; thus Msh ^ Mfof by definition. 

In FM08 and FM09 we used Mfof, the total mass of 
the particles that belong to an FOF halo by the group finder. 



Genel et al. (20081 have shown, however, that for the sub- 



set of haloes that are about to undergo minor mergers, the 
FOF mass of the lower-mass halo can be significantly higher 
(up to a factor of 1.4) than its subhalo mass {Msh) due to 
deficiencies in the FOF halo-finding algorithm. | Genel et al.| 
( [2008[ ) showed that this effect is more severe for more minor 
mergers. In light of this study, we will present the majority 
of our results for both Mfof and Msh and comment on 
the difference. 

Since the most massive haloes at 2 = are more massive 
than the most massive haloes at higher redshifts, we will 
present some results below using mass bins defined by a 
fixed percentage of haloes rather than by a fixed absolute 
mass at different redshifts. Table [l] lists the corresponding 
halo masses at 2: = 0, 0.51, 1.08, and 2.07 for each of the 
five percentile bins. Note that even in the highest 1% mass 
bin, there are ~ 5000 cluster-mass haloes available for our 
study. 



3 DIFFUSE ACCRETION 

3.1 Quantifying Halo Growth via Mergers vs 
Diffuse Accretion 



The merger trees of FOF haloes constructed in IFakhoun] 
& Ma ( 2008 1 provide, for a descendant halo of mass Mo 
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identified at redshift 20, a list of tlie masses of its Np pro- 
genitors at any cliosen zi. We label the progenitors Mi, with 
i = 1, Np, and M\ ^ M2 ^ M3... etc. The number of pro- 
genitors, Np, for a given descendant can range from 1 (i.e. no 
mergers in that timestep), 2 (binary mergers), to high values 
for massive haloes. The variable ^ is used to denote the mass 
ratio of a merging progenitor and the largest progenitor M\ : 
^ = Mi>i/Mi. 

Typically we find that the mass Mo of a descendant 
halo is not equal to the mass in the progenitors X]i>i ^'^i^ 
leading us to define 



AM = Mo 



(1) 



where AM quantifies the portion of the mass change that 
cannot be attributed to mergers of resolved haloes Mi , and 
AAf can be negative. A non-zero AM can be due to the 
merging of unresolved progenitors below our minimum halo 
mass, accretion of dark matter that is not locked up in 
haloes, or mass loss processes such as tidal stripping. We 
collectively refer to AM as diffuse accretion throughout this 
paper. 

To quantify the relative contributions to halo growth 
due to mergers vs diffuse accretion we define two mass 
growth rates for each halo of mass Mo at some redshift zo: 



Mn,cr 
Mdif 

Mtot 



At ' 
AM 
= ~Af' 

= Mmcr + Afdif = 



(2) 



Mo - Ml 
At ' 



where Mmcr is due to mergers and Mdif is due to the accre- 
tion of diffuse material. Since Mo — Mi is the mass differ- 
ence along the thickest branch of the merger tree, the sum, 
Mtot ~ Mmcr + Mdif, is simply the net mass growth rate of 
the halo. A detailed study of Aftot of the haloes in the Millen- 
nium simulation is given in McBride et al.| ( |2009[ ); here our 
focus is on the separate contributions and their environmen- 
tal dependence. The timestep At in equation (|2| is the time 
interval corresponding to Az = zi — zq, the redshift spacing 
between two Millennium outputs used to calculate the rate. 
We use Az = 0.06 for 20 = (At = 0.83 Myrs) and 0.51, 
A2 = 0.09 for 20 = 1.08, and A2 = 0.17 for 20 = 2.07, fol- 
lowing the extensive convergence tests reported in |Fakhouri| 
& Ma|((2008l). 



We also find it convenient to present the results for the 
mass growth rates in dimensionless units. When these oc- 
casions arise, we will use the fractional mass gain per unit 
redshift, dMmcr/dz/Mo and dMdif / dz / Mo, which are simply 
equal to the rates in equation Q multiplied by At/ Az/Mq. 

It is important to keep in mind that the distinction 
between mass growth from resolvable haloes and diffuse ac- 
cretion depends on the mass threshold used to define these 
two components. In this study, we have made the con- 
servative choice of including only descendant haloes with 
more than 1000 particles (with a corresponding mass of 
1.2 X lO^^M©) and progenitor haloes with more than 40 par- 
ticles (4.8x lO^^M©). A mass ratio of ^min = 0.04 is therefore 
the resolution limit on progenitors for the most poorly re- 
solved descendant haloes (1000 particles) in our sample. For 



these haloes, the mass contributed by progenitors above 40 
particles is tagged as due to mergers, while any remaining 
mass contribution is tagged as due to diffuse accretion. 

For more massive descendants, there are two natural 
choices for the mass threshold: one can include either all pro- 
genitors down to the 40 particle resolution limit as mergers, 
or all progenitors down to a fixed mass ratio of ^min = 0.04 
as mergers (any progenitors with ^ < ^min are then counted 
as "diffuse" material). The former has the advantage that 
40 particles correspond to a smaller ^min for more massive 
descendants; we will therefore have better statistics and dy- 
namic range for more massive haloes. The latter choice has 
the advantage that mergers and diffuse accretion can be 
compared objectively across different descendant mass bins 
since the two components are defined with respect to the 
same ^min for all haloes and are, therefore, effectively equally 
well resolved. In Sec. 15.11 we will show that the basic be- 
haviour of the environmental dependence presented in this 
paper does not depend on which definition is used. We gen- 
erally favour the latter definition, and will use ^min = 0.04 
to separate mergers vs. diffuse accretion unless otherwise 
stated in this paper. 

3.2 Distributions of Afmcr and A/dii 

Fig.[l]presents the distributions of the mass growth rates due 
to mergers (left panel) and diffuse accretion (right panel) at 
2 = for five halo mass bins. The two halo mass definitions 
discussed in Sec. 2.3, Mfof (dashed) and Msh (solid), are 
compared. In constructing these plots, we have used a ^ > 
^min cut: mergers with mass ratios less than S^rnia ~ 0.04 are 
counted as diffuse accretion. 

We have also imposed a cut on the very negative and 
positive tails of the distributions. Although the majority 
of the haloes in the Millennium simulation have reasonable 
dMmci/dz/Mo and dMdn/dz/Mo, about 0.3% of them more 
than doubled their mass via diffuse accretion between 2 = 
and 2 = 0.06, and about 0.4% lost over half of their cur- 
rent mass via diffuse accretion. Our detailed examination of 
the merger histories of haloes in these tails of the distribu- 
tions shows that these very large values of |dMdif/dz/Mo| 
are unphysical and are due primarily to complex and spu- 
rious halo fragmentation. The algorithm used to construct 
the FOF merger trees, stitch-3, only stitches halo fragmenta- 
tions that are destined to remerge within 3 snapshots. Some 
fragmentations do not satisfy this condition. 

This small set of haloes - a set too small to affect 
the merger statistics studied in our earlier work - can 
significantly affect the mean values of dMmcr/dz/Mo and 
dMdii/dz/Alo- This is especially true when these mean 
values are computed as a function of halo environment, 
since fragmentations occur more frequently in the messier 
dense environments. Moreover, haloes with artificially high 
growth rates have correspondingly artificially low forma- 
tion redshifts (2/). These events pollute the distribution of 
2/ and tend to underpredict the formation redshift, espe- 
cially in overdense regions. We therefore implement a cut 
on dMuiei/dz/Mo and dM^a/dz/Mo to remove this aber- 
rant population. This is done by demanding that a halo 
never gains more than a fraction /+ = 1/2, nor lose more 
than a fraction /_ = 1/2, of its final mass Mo in either 
dMmer/dz/Mo OT dMdi! / dz / Mo between any two adjacent 
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Figure 1. Distributions of the mass growth rates (per unit redshift) via mergers (left panel) and via diffuse accretion (right panel) for 
the ~ 500, 000 haloes at 2: = in the Millennium simulation. In each panel, the colours denote haloes of different mass percentiles (see 
Tablejljl, and two definitions of halo mass are shown (dashed for MpoF\ solid for Mgjj). 



simulation outputs along the lialo's main branch (dotted 
vertical lines in Fig. [T]). This is a fairly stringent cut as it 
is apphed to a halo's entire mass history. The resulting cut 
removes ~ 3% of the haloes in our sample. Varying the val- 
ues of /+ and /_ in the range of / € [0.2, 1.] and allowing 
/+ 7^ /_ made insignificant changes and did not affect the 
overall trends reported in this work. A handful haloes with 
rates beyond these cuts remain in Fig. [l] because the cut is 
applied to the halo growth across adjacent outputs whereas 
the rates shown in Fig. [l] are computed across the three 
outputs between z = and z — 0.06 (which straddle the 
z = 0.02,0.04 outputs). 



In the left panel of Fig. [T] the strong peak at 
dMmcr/dz/MQ = is due to the fact that between z = Q 
and 0.06 roughly 65% of haloes have only one progeni- 
tor and, therefore, experience no mergers. The remaining 
35% of haloes have a distribution that falls off smoothly 
out to dMraci / dz / Mq ~ 8.3 (dotted vertical line). Values 
exceeding 8.3 correspond to haloes that have gained more 
than half their mass between z — 0.06 and z — (i.e., 
Mi/Mo/ Az > 1/2/0.06 = 8.3); only ~ 0.2% of haloes 
are in this category. In the right panel of Fig. [l] the peak 
occurs at small positive values of dMdn/dz/Mo, though a 
significant set of haloes experience diffuse stripping. 



To complement the differential distributions presented 
in Fig.[l] we show in Fig. [2] the cumulative distributions for 
the number density of haloes above a given mass growth 
rates in units of Mq per year. The three panels are for 
three halo mass bins. The contributions from mergers (thick 
curves), non-halo material (thin curve), and the sum (black 
curves) are plotted separately. The curves for the total rate 
are identical to the solid curves (z = 0) in the bottom panels 
of Fig. 5 oflMcBride et al.l (120091). This figure shows that the 



haloes with high accretion rates in the simulation are under- 
going mergers with other haloes, while the haloes with low 
accretion rates are mostly accreting non-halo mass. 



3.3 Redshift Evolution of the Mean Growth Rates 

Having examined the distributions of the growth rates a,t z = 
in the last section, we now study the redshift dependence 
of the mean growth rates. The left panel of Fig. [3] shows 
the mean Mmor (solid), Mdif (dashed), and total growth rate 
Afmcr + A^dif (dotted) as a function of redshift for three mass 
bins: lO^^M© (blue), lO^M© (green), and IO^A^q (red), 
where the threshold ^min = 0.04 is used to define Mmor and 
Mdif. The right panel shows the same information except the 
rates are expressed in the dimensionless units of dM/dz/AIo- 
The growth rates per year (left) are seen to increase rapidly 
with increasing z, while the rate per unit z (right) has a very 
weak dependence on z out to z ~ 6. As a function of halo 
mass Mo, the right panel shows that the haloes of higher 
mass experience somewhat higher fractional mass growth 
rates than haloes of lower mass. Similar dependence on time 
and mass was seen for the halo merger rates in |Fakhouri &:| 
MalpOOSl. 



The left panel of Fig. [3] can be compared directly with 
the total mass growth rates shown in Fig. 4 of |McBride et al.] 
( 2009 1, which is well approximated by equation (7) proposed 
there: 



Mtot\ = 42MQyr-' 



M 



\ 1.127 
IOI2M0J 

x(i + i.i7z)^nm{i + zy + 



(3) 



where Qm and Qa are the present-day density parameters 
in matter and the cosmological constant, and we have as- 
sumed Qrn + = 1 (used in the Millennium simulation). 
We have overlaid this fit using black dotted curves in Fig. [3] 
A comparison of the dashed and solid curves in Fig. [S] shows 
that the mean Mdif is consistently ~ 30% higher than the 
mean A^nicr at all mass and redshifts. The functional form of 
the mass and redshift dependence in equation ([3| is there- 
fore also applicable for these separate rates; only the overall 
amplitude needs to be adjusted to obtain a fitting form for 
Mnicr and Mdif. 

The similarity in shape and mass dependence of the 
Afmor and Afdit curves in Fig. [3] may lead one to conclude 
that the "non-halo" material contributing to Mdif can be 
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dM/dt [Mq /yr] dM/dt [Mq /yr] dM/dt [Mq /yr] 

Figure 2. Cumulative distributions of the z = mass growth rates via mergers with haloes (thick curves) and diffuse accretion (thin 
curves) for haloes in three mass bins (left to right panels). The axes are in physical units: halo number density (y-axis) and accretion 
rate in Mq per year (x-axis). The total rate is plotted in black. The solid and dotted curves are for the two halo mass definitions Mgu 
and MpoF- The high accretion rate events are dominated by mergers with other haloes. 
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Figure 3. Redshift evolution of the mean halo mass growth rates due to mergers (solid) and diffuse accretion (dashed) for three mass 
bins: IO-'^^Mq (blue), IO^^Mq (green), and 10"*^ ''Mq (red). The rates in units of Mq per year (left panel) are seen to rise with increasing 
2, while the rates in units of per redshift (right panel) is largely independent of z. We also plot the total merger rate (haloes + diffues) 
with dotted lines and overlay the fit presented in equation[3]in black. Note that the fit is quite good and the total rate and fit are difficult 
to distinguish. 



simply attributed to mergers with sub-resolution haloes. 
This interpretation is too simplistic, however, as we will dis- 
cuss in Sections 4 and 5 below where their environmental 
dependence is analysed. 

4 THE ENVIRONMENTAL DEPENDENCE OF 
HALO GROWTH RATES AND HISTORIES 

4.1 Halo Mass Growth Rate due to Mergers vs 
Diffuse Accretion 

In Part I of this series ( jFakhouri fc Ma]|2009| |, we reported 
a strong positive correlation between the number of halo 



mergers and halo environment. We now turn to the environ- 
mental dependence of the rate of halo mass growth. Since 
haloes acquire mass through both mergers and diffuse ac- 
cretion in the form of non-halo mass or unresolved haloes, 
we examine separately the two growth components defined 
in equation ([2|. 

The mininnim mass ratio ^min — 0.04 discussed in p?T] 
is used as the threshold between Mmer and Mdn- We com- 
pute Mmci and Mdif for each descendant halo and then bin 
the results by averaging over various ranges of mass and 5. 
For haloes with only one progenitor (above the mass thresh- 
old), there is no contribution from mergers by definition, and 
Mmer = 0. These events are included when averaging over 
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Figure 4. Environmental dependence of the halo mass growth rates due to contributions from merging progenitor haloes (top), diffuse 
accretion of non-halo mass (middle), and sum of the two (bottom). Four redshifts z = 0, 0.51, 1.08, 2.07 are shown for comparison (left to 
right). Within each panel, the five curves show five mass bins for the descendant haloes listed in Table[l](blue: low mass, red: high mass) 
with solid lines computed using the subhalo mass definition and dashed lines computed using the FOF mass definition. On average, 
haloes in denser regions experience a higher mass growth rate from mergers and a lower rate from diffuse accretion than similar mass 
haloes in the voids. These two opposite trends with environment roughly cancel out to yield a weak, but net negative, dependence on 
1 + 57_FOF for the total growth rate. 



Mmcr SO that a fair comparison can be made to the mass 
growth rate due to diffuse accretion, Mdu- 

Fig. |4] compares the environmental dependence of the 
mean halo mass growth rate (per unit redshift) due to merg- 
ers (top row), diffuse accretion (middle), and the total rate 
(bottom) computed from the Millennium simulation. The 
rates at four redshifts are shown: z — 0,0.51, 1.08 and 2.07 
(from left to right). Within each panel, different colours cor- 
respond to different mass bins (listed in TablejlJ , where more 
massive haloes have higher growth rates. The solid curves 
are computed using AIsh , the sum of all the subhalo masses 



in a given FOF halo, while the dotted curves are computed 
using Mfof (see Sec. 2.3). Both mass definitions yield sim- 
ilar quantitative behaviour, though the subhalo mass defini- 
tion yields slightly lower merger growth rates in accordance 



with the results in Genel et al. ( 2008 1 



The mass growth due to mergers (top row) is seen to 
increase in denser regions, in agreement with the higher 
merger rates in denser regions reported in FM09. This is 
to be expected as the merger rate B/n is related to the 
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mass growth due to mergers Mmer by 

nl 



0.04 



(4) 



per- 



As with B/n, the environmental dependence of AI„ 
sists out to higher redshift. 

By contrast, the middle row of Fig. |4] shows that the 
mass growth rate due to non-halo material has an opposite 
dependence on 5: the mean value of dMm/ dz / Mq decreases 
with increasing 5 for all halo masses at all redshifts. This 
trend is clean when Msh is used as the mass definition (solid 
curves) but the curves show a sharp rise for low mass haloes 
when Mfof is used (blue dotted curve). We believe this 
is due to the inadequacy of the FOF halo finder in dense 
regions. As Genel et al. ( 20081) have shown, the FOF mass for 
low mass haloes sometimes rise as they approach high mass 
haloes; this may contribute to the rise in FOF mass observed 
in the densest regions. Despite this fact, we emphasise that 
the discrepancies between the Mfof rates and the Msh 
rates only arise in the densest regions (beyond 1 + (57_fof ~ 
5). Outside of these regions, dM^ci/dz/AIo still increases 
with 5 and dMdn / dz /Mo decreases with 5 regardless of the 
halo mass definition. 

The mean total mass growth rate is shown in the bot- 
tom row of Fig. |4] It shows that the strong but oppo- 
site dependence on environment found for the two growth 
rates act against each other, producing a relatively weak 5- 
dependence for the overall growth rate that declines with 
increasing 5 for the three low mass bins. This dependence 
is not as clearly seen for the very highest mass bin, though 
the statistics are poor for the cluster-scale bin (red curves). 
Again, the trends are much cleaner for Msh than for Mfof 
in the high-(5 regions. The mass dependence of the total 
growth rates, on the other hand, is more prominent than 
the environmental dependence, with massive haloes grow- 
ing more rapidly than low mass haloes. Equation (7) of 
McBride et al. ( 2009 \ (repeated above in eq. ([3| shows that 
this mass dependence is well approximated by a power-law: 

The curves in Fig. |4] are sensitive to the choice of cut 
used to remove the aberrant haloes (see Sec. 3.2) at the 20- 
50% level. The qualitative features of the plot, including the 
markedly opposite correlation of the two rates with 5, and 
the negative, but relatively weak, S dependence of the total 
rate, are insensitive to choices of /+ and /_ in the cut and 
persist even if the cut is made asymmetric (/+ 7^ /_). 

It should be pointed out that although the two growth 
rates depend on the choice of ^mi„ used to define resolv- 
able haloes and diffuse material (0.04 in this case), the total 
growth rate is independent of this parameter. Thus, redis- 
tributing progenitors with ^min < 0.04 from dAddn/dz/Mo 
to dMmai/dz/Mo would not impact the bottom row of Fig.|4] 
We will return to this point in more detail in §5.1| 



4.2 



Fraction of Final Halo Mass Acquired from 
Mergers vs Diffuse Accretion 



In addition to the instantaneous halo mass growth rates 
shown in Fig. |4j a related quantity of interest is the inte- 
grated contribution to a halo's final mass due to mergers vs 
diffuse accretion. Fig. [5] compares these contributions as a 
function of environment. The lowest set of curves plots the 
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Figure 5. Environmental dependence of the fraction of the final 
halo mass Mq gained by mergers vs diffuse accretion. The bottom 
two sets of curves show the rising mass fraction with 1 -|- S^—fof 
for mass gained via mergers with progenitors of mass ratios ^ > 
1/3 and ^ > 0.1. The top two sets of curves show the declining 
mass fraction with 1 -|- ^y-poF for the contributions from non- 
halo mass with ^ < 0.04 and g < 0.1. Within each set of curves, 
different colors represent different mass bins, and solid and dotted 
curves are for the two halo mass definitions Msh ^'Ud MfoF} 
respectively. This figure illustrates the increasing importance of 
mergers to halo growth in overdense regions and importance of 
diffuse accretion in the voids. 



fraction of a halo's final mass gained from major mergers 
(defined to have progenitor mass ratio ^ > 1/3 at the time 
of merger) through the halo's entire history (out to z = 6); 
the set of curves above it shows the same quantity for more 
minor mergers with mass ratio ^ > 1/10. The five curves 
in each set show the five halo mass bins listed in Table [T] 
Similar to the top panels of Fig. |4j these two sets of curves 
exhibit a positive environmental dependence, where haloes 
in denser regions experience more mergers and have a higher 
fraction of their mass coming in from mergers. From voids to 
overdense regions, the average mass fraction of a halo due to 
major mergers (^ > 1/3) increases from ~ 20% to 25%, and 
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the mass fraction due to mergers with ^ > 1/10 increases 
from ~ 30% to 40%. 

The upper two sets of curves show the analogous quan- 
tity for the fraction of a halo's final mass gained from ma- 
terial with ^ < 0.04 and ^ < 0.1. In contrast to the lower 
two sets of curves, the correlation of the mass contribution 
by non-halo matter with 5 has the opposite sign, again con- 
sistent with the middle panels of Fig. |4] 

A consistent picture is therefore emerging in that haloes 
in denser regions experience a higher rate of growth via 
mergers and a lower rate of growth via diffuse accretion than 
haloes of comparable mass in less dense regions. As a result, 
in dense regions, a higher fraction of a halo's final mass is 
gained through mergers than haloes in the voids. By con- 
trast, diffuse accretion contributes less to halo mass growth 
in denser regions than in the voids. Overall, however, diffuse 
accretion (defined as ^ < 0.04) is an important component 
of a halo's final mass regardless of the environment: its con- 
tribution to the mass fraction is never lower than ~ 40% 
(second set of curves from the top in Fig. [5|. 

4.3 Halo Formation Redshifts 

We follow the convention in the literature and define the for- 
mation redshift Zf for a halo to be the redshift at which the 
halo's mass first reaches Mo/2 (when traced backwards in 
time), where Mq is measured at z = 0. This is computed by 
tracing the mass of each of the ~ 500, 000 haloes from low 
to high redshift in the Millennium outputs. A complemen- 
tary Zf can be defined by going from high to low redshift; 
we have done this as well and found little difference in our 
results. 

Fig. |6] shows the mean formation redshift for haloes in 
five mass bins (in increasing mass from top down) as a func- 
tion of their environment at 2; = 0. The two panels compare 
two definitions of the local densities discussed in IFakhounl 
|fc Ma| ( [2009[ ) : I + S^y-fof used throughout this paper (right) 
and the simpler 1-1-^7 (left). The two measures of local 
densities give very similar results except in the most mas- 
sive halo bin (red curves), where Zf shows little correlation 
with ^T-FOF but decreases with increasing 87. The negative 
correlation of Zf with 67 should not be interpreted as a true 
dependence on a halo's surrounding environment; instead, it 
is due to the tight correlation between 67 and the masses of 
massive haloes shown in Fig. 2 of FM09 and the fact that, on 
average, more massive haloes form earlier than less massive 
haloes. This degeneracy between mass and density of mas- 
sive haloes is removed when the variable 57_fop is used (see 
Fig. 2 of FM09). The resulting Zf vs Jt-fof shows almost 
no correlation in the right panel of Fig. [6] 

The solid and dotted curves in Fig. jGlcompare the mean 
values of 2:/ computed using the Msh vs Mfof definition of 
halo mass. As seen in Fig. |4j the two mass definitions yield 
similar results in underdense regions but differ in the dens- 
est regions. Fig. |6] also compares Zf computed with (darker 
shades) and without (lighter shades) the cut described in 
Section |3.2[ The two z/'s are identical in underdense re- 
gions, but the sample without the cut has a lower Zf aX high 
5. This drop is unphysical and is due to the orhpaned halo 
fragments that suddenly appear at low z and thus bias Zf 
to artificially low values. 

A number of conclusions can be drawn from Fig. |6] 



First, within a given environment, higher mass haloes on 
average form later (i.e. lower zf). For instance, the mean for- 
mation redshift drops from 2/ ~ 1.1 to 1.4 for M ~ IQ^^Mq 
haloes, to about 0.7 for M > W^Mq haloes. This is a weU 
known property of hierarchical cosmological models such as 
the ACDM. This mass dependence is also consistent with 
that shown in the bottom panels of Fig. |4] where more mas- 
sive haloes have a larger mean mass growth rate and hence 
a later mean formation time. 

The second feature to note is that galaxy-size haloes 
(top few curves) in denser regions form earlier than the 
same mass haloes in the voids, whereas zj for cluster-size 
haloes (bottom curve) hardly correlates with environment 
(when 57_FOF is used as a measure). The earlier formation 
of galaxy haloes in denser regions is consistent with the halo 



assembly bias discussed in a number of recent papers (|Gao 
ct al. 2005; Shcth & Tormcn 2004 Jing et al. 2007; Gao & 


White,2007,, Barker et al. 


2006,,Wechsler et al.,2006,, Wetzel 


et al. 


2007 


|Wang et al. 


20071 |Hahn et al 


|2008 McBride 


et al. 


2009 


1. Interestingl;y 
1 reported deer 


Wechsler et al. 


2006 


1 and Jing 


et al. 


2007 


easing zj with increasing halo bias 



for 10 



10^ 



Mq haloes. This trend is similar to the trend 
for the cluster-mass haloes (red curve) in the left panel of 
our Fig. [6] Since the connection between 5 and halo bias is 
complex, it is not immediately clear that their results are 
inconsistent with ours. Perhaps the negative correlation ob- 
served by these authors is also due to the fact that the bias 
of high mass haloes is strongly correlated with mass and 
does not provide an independent measure of environment. 

The third feature to note in Fig. |6] is the slight upturn 
of 2/ in the most underdense void region. This trend makes 
intuitive sense since haloes in sufficiently underdense regions 
are deprived of fuel for growth, thereby growing more slowly 
and having larger formation redshifts. This trend does not 
seem to appear in any of the earlier literature exce pt |Des- 
|jacques| ( [2008| ), which noticed a similar trend in ^Harker et al. 
( 2006 1 . We suspect that the use of bias as a measure of halo 



environment in most of these earlier papers (as opposed to 
the local density used here and in Harkcr et al. 2006) is not 
a sensitive probe of the 2/ statistics of haloes in the voids. 



4.4 Mass Reservoir outside Haloes 

In this subsection, we address the question: are the opposite 
environmental correlations of Afmcr and Mdif with 5 seen in 
Figs. [4] and [5] also present in the surrounding mass reservoirs 
outside of the haloes? If yes, then the different environmental 
dependence of the growth rates is originated from the mass 
reservoirs that are feeding the haloes. If not, then different 
accretion timescales have to be operating for the haloes and 
non-halo material in overdense vs underdense regions. The 
answer to this question will therefore help explain the origin 
of the difference in the environmental correlations reported 
thus far. 

To answer this question, we analyse the halo vs non- 
halo composition of the mass in the volume used to compute 
57_FOF for each halo at 2 = 0. Centered at each halo, we 
compute the mass within the 7 Mpc sphere (excluding the 
FOF mass itself) that resides in haloes with mass ratios 
exceeding ^ > 0.04 (relative to the central FOF mass) vs 
in all other matter. The results are shown in Fig. [7] which 
plots the mean values of these two components in the mass 
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Figure 6. Environmental dependence of the mean formation redshifts for haloes in five mass percentile bins (same bins as Fig. [sjl in 
increasing mass from top down. Two measures of environment are shown for comparison: I + S7 that includes the halo mass (left), and 
1 + Sj_pQp that excludes the central halo. Solid lines are computed using Mgu and dotted lines using MpoF- The darker colors are 
computed after the cut described in Section |3.2| is applied; the light colors are computed without the cut. This figure shows that more 
massive haloes (bottom curves) on average form more recently. Within a mass bin, galaxy-size haloes in denser regions form earlier than 
those in the voids, while cluster-size haloes show little environmental dependence in their formation time. 
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Figure 7. Composition of the mass reservoir in the region outside 
of the halo but within a 7 Mpc sphere centered at the halo. The 
vertical axis shows the contribution to l + (57_poF by haloes (with 
mass above 4% of the central halo; solid curves) and the remaining 
material (dashed lines). By construction, the solid and dashed 
curves add up to the thick diagonal line. Different colours denote 
the same mass bins in earlier figures. The mass reservoir for haloes 
in overdense regions is seen to be dominated by resolved haloes, 
whereas the reservoir for haloes in the voids is dominated by non- 
halo material. 



reservoir as a function of 1 -I- 57_fof for five mass bins. The 
black diagonal line marks the sum of the two components, 
which by definition is simply 1 + St-fof- 

Both the halo and non-halo components in the 7 Mpc 



reservoir outside the haloes are seen to increase with 1 -|- 
57_FOF, but the halo component has a much steeper slope. 
The two components contribute a comparable amount for 
haloes in the cosmic mean density (i.e. 57_fof ~ 0), but up 
to 80% of the mass reservoir is made of haloes in very dense 
regions, while this fraction drops to ~ 20% in voids, where 
the diffuse component dominates. Even though only a frac- 
tion of this mass reservoir is accreted onto the central haloes 
and contributes to actual halo growths, Fig. [7] suggests that 
the higher merger rates in denser environments are a direct 
consequence of the larger fraction of available haloes in the 
reservoir. 

A somewhat subtle but important feature to note in 
Fig. [7] is that the absolute amount of diffuse mass in the 
reservoir rises with increasing 57_fof even though the frac- 
tional contribution due to this component drops. This rising 
trend with 57_fof of the diffuse component is in striking 
contrast to the middle panels of Fig. [4] where the accre- 
tion rates of diffuse material, Mdif, onto haloes are seen 
to decrease in denser regions. The naive interpretation that 
Afdif decreases with increasing 57_fof is due to a dwindling 
reservoir of diffuse material outside of the haloes in denser 
regions is therefore invalid. On the contrary, at a fixed halo 
mass, there is more non-halo material available in the re- 
gions where Mdif is lower. 

Additional gravitational effects are likely to be at work 
in dense environments to explain why the halo mass growth 
rate due to diffuse material is stunted in denser regions while 
the supply of this mass is in fact larger. For instance, the 
non-halo component in the reservoir may be dynamically 
hotter than the resolved haloes due to tidal heating and 
stripping, thereby reducing Mdif in higher density regions as 
shown in Fig. |4] Pieces of this picture have been suggested 
in the literature. [Wang et al. (20071 noted that a halo's 
ability to consume the available fuel is as important as the 
amount of fuel available for growth. They found that haloes 
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in dense regions accreted matter more slowly than expected, 
although they did not examine the separate A/mor and Mdn 
components as we do here. It has also been suggested that 
low-mass haloes form earlier because their late accretion is 



eliminated in a competition with massive neighbors (Zent- 
ner|2007 |. Fig.[4j however, shows that this competition idea 
is not strictly true: regardless of mass, all haloes in over- 
dense regions have both increased growth rate Mmor and 
decreased growth rate Mdu- The growth rates of low mass 
haloes are therefore not simply controlled by competition 
with nearby high-mass haloes in overdense regions. Alterna- 
tively, the diffuse material surrounding very dense regions 
may preferentially accrete onto the lower-density filaments 
and the haloes within them, which then infall onto the clus- 
ters. Further analysis taking into account of filaments as well 
haloes would be required to test this idea. 



4.5 Time Evolution of a Halo's Environment 

Thus far we have used the mass densities computed at a 
halo's redshift to quantify its environment. Here we exam- 
ine whether the environments of haloes evolve significantly 
over their lifetimes, that is, if the haloes in overdense re- 
gions today tended to have progenitors that also resided in 
overdense regions at earlier times. 

To study such environmental evolution, we first group 
the haloes a,t z — into mass bins. Within each mass bin, 
the haloes are further divided into five environmental bins, 
each corresponding to 20 percentile in the distribution of 
1 + (Jt-fof. Each halo is therefore assigned to a mass and a 
57_FOF bin. For each halo at 2: = we then identify its most 
massive progenitor at z = 0.51 and z — 1.08. We assign 
each progenitor a mass bin at its redshift and compute the 
(Jy-FOF bin to which it belongs, using the density field at 
that given z. 

We then compute the fraction of haloes at z; = 0, as 
a function of mass bin, whose most massive progenitors re- 
side in the same percentile S bin at z = 0.51 and 1.08. Such 
haloes have not deviated from their local environment sig- 
nificantly over time. Since 5 evolves with time, it is more 
meaningful to compare relative percentiles instead of abso- 
lute values of S at different redshifts. We also compute the 
fraction of haloes at 2: = that are at most one percentile 
bin away a,t z — 0.51 and 1.08 . These haloes have moved 
slightly outside of their z = environmental percentile. 

We find that ~ 70% of M < 10"Mq haloes do not 
change environmental bins between z = 0.51 and today, and 
60% do not change environments between z = 1.08 and to- 
day. In addition, more than 95% of these haloes are within 
one 5-percentile bin of their final environment. The range 
of S for cluster-mass haloes is too narrow (see, e.g.. Fig. 6) 
for this analysis to be meaningful. Thus, the majority of 
haloes reside in the same environmental context across their 
lifetimes. 



and diffuse non-halo material (Afdif), and their respective 
environmental dependence. As emphasised there, merging 
progenitors with a mass ratio of ^ > ^min = 0.04 were 
counted towards Mma , whereas the rest of the mass growth 
was counted towards Mdu- 

To test the robustness of the results in Sections 3 and 
4 with respect to the value of ^min used to define haloes 
vs non-haloes, we take haloes in one of the more massive 
bins (the 90-99% mass bin, corresponding to 1.4 x 10^'^ to 
1.1 X lO^^AfQ) that provides good statistics as well as high 
mass dynamic range, and recompute the two rates using 
different values of ^min. The results are shown in Fig. |8j 
which plots the environmental dependence of the two mass 
growth rates for three values of ^min: 0.4 (black), 0.04 (dark 
grey; the value used in earlier figures), and 0.004 (light grey). 
As ^min is lowered, more of the mass is counted towards the 
merger component, so dMmci/dz/Mo (solid curves) increases 
and dMdii/dz/Mo (dashed curves) decreases, while the sum 
of the two (dotted curve) is, by construction, independent 
of ^min. We note that even though the overall amplitudes 
of the two rates change with ^min, dM^ncr / dz / Mq remains 
positively correlated and dMdu / dz/ Mq remains negatively 
correlated with 57_fof regardless of ^min. This test suggests 
that the environmental correlations shown in Figs. |4] and [5] 
are robust to the value of ^min used to define haloes vs non- 
haloes. 

We explore further the nature of the diffuse non-halo 
component by examining the hypothesis that this compo- 
nent is comprised entirely of sub-resolution haloes such that 
in the limit of ^min (or more precisely, the true mini- 
mal halo mass, which is likely to be about an Earth mass 



(Diemand et al. 20051, but the difference is negligible), we 
would expect Afdif — > and Mma —* Mo — Mi. To assess 
whether this hypothesis is plausible, we recompute Mmei by 
extrapolating down to ^min = 0, assuming that Mma main- 
tains the same power-law dependence on ^ determined for a 
given mass bin in the simulationj^ 

The predicted mass growth rate under this assumption 
is plotted as a function of environment in Fig. [s] (top solid red 
curve). In comparison with the total rate Mma + Mdn (dot- 
ted curve) determined from the Millennium simulation, we 
observe a clear environment-dependent gap between these 
two curves that widens in the lower-density regions. This gap 
suggests two possible scenarios: (1) the power-law behaviour 
of Mmor steepens towards low ^ (an effect that is not ob- 
served down to ^ ~ 10~* resolvable by current simulations), 
or (2) there exists a truly diffuse component of dark mat- 
ter that does not belong to bound haloes. Moreover, both 
scenarios must be dependent on environment. For (1), the 5- 
dependent gap in Fig. [S] can only be closed if the power-law 
slope of Mmar VS ^ for sub-resolution progenitors is a strong 
function of halo environment (being steeper in underdense 
regions). We cannot rule out this possibility without testing 
it in higher-resolution simulations, but all the halo proper- 
ties that we have examined thus far in FM08 and FM09 



5 DISCUSSION 

5.1 Interpreting "Diffuse" Accretion 

In Sections 3 and 4, we presented results for the mass ac- 
cretion rates of dark matter haloes through mergers (Afmcr) 



^ We note that our best-fit slope for the ^-dependence of Afmer 
in [Fakhouri &: Ma| l |2008| was in fact slightly steeper than —2, 
which would give a divergent growth rate. Our goal there was to 
obtain a simple universal form. For better accuracy, here we fit 
A^mor vs g for each mass bin independently, where the slope had 
a very mild mass dependence, ranging from —1.7 to —1.9. 
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Figure 8. Dependence of the two growth rates dMmcr/dz/Mo 
(solid) and dM^if / dz/ Mq (dashed) on the mass threshold ^min 
used to define the halo and non-halo components. The sum of 
the two rates (dotted), by construction, is independent of ^niin- 
Haloes in the 90 — 99% mass bin are used. As ^min is lowered from 
0.4 (black), 0.04 (dark grey), to 0.004 (light grey), the amplitude 
of dM^ar/d,z/Mo increases and that of dM^if/ dz/ Mq decreases, 
but their respective correlations with <57_pop remain the same, 
suggesting that the environmental dependences reported in this 
paper are robust and do not depend on the exact values of ^min 
used to define haloes vs non-haloes. The top solid curve (in red) 
shows the hypothetical value of dMmer/dz/Mo extrapolated down 
to 5min = (see text for discussion). One explanation for the 5- 
dependent gap between this curve and the total rate (dotted) is 
that there exists a truly diffuse component of dark matter that 
does not belong to bound haloes, in particular in the lower-density 
regions. 



do not exhibit such complex, non-universal behavior. On 
the contrary, down to the current simulation resolution, the 
halo-halo merger rates and mass growth rates as a function 
of progenitor mass ratio ^ were all well fit by a single power 
law, where the slope was entirely independent of the halo 
environment (see Figs. 4, 5, and 7 of FM09). We therefore 
favour scenario (2) that the diffuse non-halo mass has a truly 
diffuse component, and this component is more prominent 
in lower-density regions. New simulations with higher mass 



resolution such as the Millennium-II ( Boylan-Kolchin et al. 



2009 1 will be used to explore further these two possibilities. 



5.2 



Implications for the Extended 
Press-Schechter and Excursion Set Models 



The Press-Schechter model (Press & Schechter 19741 as- 



sumes that primordial Gaussian-distributed dark matter 
perturbations collapse into haloes of mass M when their 
linearly extrapolated overdensities, smoothed on a scale cor- 
responding to a given mass, exceed a critical threshold com- 
puted using a simple spherical collapse model. 

This framework was extended to address the problem 
of halo growth in the excursion set formalism (sometimes 



called the Extended Press-Schechter or EPS model) (Bond 
et al. 1991 Lacey & Cole 19931. Using a fc-space tophat 



smoothing window, these authors showed that the density 
perturbation at a given point in space undergoes a random 
walk as the perturbation is smoothed on smaller and smaller 
scales. These perturbation "trajectories" obey the diffusion 
equation and are Markovian, i.e., the change in overdensity 
at a given scale is independent of the overdensities at other 
scales. The EPS model in its standard form therefore pre- 
dicts no environmental dependence. For example, the con- 
ditional mass function (j>{Mi, zi\Mo, zq) - the distribution 
of progenitor mass Mi at redshift zi for a descendant halo 
of mass Mo at zq - is predicted to be the same for all en- 
vironments, in contrast to Fig. 7 of Fakhouri & Ma ( 2009 1 



that showed descendant haloes in dense regions to have more 
progenitors than those in voids. 

Recent attempts have been made to introduce environ- 
mental correlations into the EPS model. As emphasised in 
Zentner (20071, the Markovian nature of the random walks 



in the excursion set model is not a prediction but is an as- 
sumption originating from the use of the fc-space tophat win- 
dow function. This window function is traditionally chosen 
to reduce the time required for computing density perturba- 
tion trajectories since for other window functions, the per- 
turbations must be computed at all scales simultaneously. 



When a Gaussian window function was used, Zentner ( 2007 \ 



indeed found an environmental dependence of the halo for- 
mation redshift, but the dependence was opposite to that 
seen in numerical simulations (including Fig. |6] of this pa- 
per), where older (i.e. earlier forming) haloes tended to be 
more clustered. 



Desjacques (20081 instead chose to introduce environ- 



mental dependence into the threshold density for collapse in 
the ellipsoidal excursion set formalism. Overdense regions 
will tend to exert more tidal shear on collapsing haloes, an 
effect that causes the haloes to virialize earlier. Like IZent-l 



(20071, however, Desjacques (20081 found the low mass 



haloes in denser regions to form later than similar mass 
haloes in emptier regions, again opposite to the effect seen 
in simulations. Desjacques (20081 did note that his results 



may be valid in the lowest-density region, where Fig. |6] here 
and Harker et al. ( 2006 1 both reported a subtle increase in 



formation redshift with decreasing 5. 



Sandvik et al. (20071 carried out a multi-dimensional 



generalisation of EPS that allowed them to incorporate envi- 
ronmental effects by tracking the shapes of collapsing haloes 
and incorporating information about gravitational shear. 
They found 2:/ in this approach to depend very weakly on 
the environment; moreover, the S dependence they did ob- 
serve was stronger for more massive haloes, in contrary to 
the results of N-body simulations such as our Fig. [6] They 
proposed, instead, that the environmental dependence of a 
halo's formation history may be related to the halo's progen- 
itor pancakes and filaments. They found that haloes whose 
progenitors were the most massive pancakes or filaments 
were more clustered than halo's whose progenitors were the 
least massive pancakes or filaments, and that this clustering 
dependence was stronger for low mass haloes. 

In summary, there is not yet a modification of the EPS 
model that can successfully predict the correct correlation of 
halo formation with environment seen in numerical simula- 
tions. We emphasize that in addition to the "assembly bias" 
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for the formation redshift discussed in Sec 4.3 and shown in 
Fig. |6j an improved EPS model would also need to predict 
the strong but opposite environmental dependence of the 
growth rates shown in Fig. [4] A first step towards this goal 
would be to introduce into the model, and track, the con- 
tribution made by the non-halo material Mdif. The model 
would also need to be sophisticated enough to account for 
the opposite behaviours of Mmer and Mdif as a function of 
the local density. 



6 CONCLUSIONS 



We have used the Millennium simulation (Springel et al 
|2005[ ) to investigate the dependence of halo mass growth 
rates and histories on halo environment. This paper comple- 



ments our previous paper (Fakhouri & Ma 20091 in which 



we reported and quantified how haloes in overdense regions 
experienced higher merger rates than haloes in underdense 
regions. Here, we have studied the mass growth rates due to 
both mergers with other resolved haloes, Mmor, and accre- 
tion from non-halo "diffuse" material, Mdif. Results for the 
distributions of Mmcr and Mdif at 2; = and the redshift 
evolution of the mean Mmer and Mdif are summarized in 
Figs.[l][3l 

As a function of halo environment, we have found (see 
Fig.|4j| that the growth rate Mmer due to mergers, in agree- 
ment with the merger rates in Fakhouri & Ma ( 2009 1 , is 



positively correlated with the local density, whereas Mdif is 
negatively correlated. Consequently, in denser regions, merg- 
ers play a relatively more important role than diffuse accre- 
tion for the mass growth of haloes, and a higher fraction 
of a halo's final mass is acquired through mergers in these 
regions than in the voids (Fig. [sj. 

We have shown that the origin of the environmental de- 
pendence of growth rate due to mergers, Mmor, is directly 
linked to the mass reservoir immediately outside the virial 
radii of the haloes (see Fig. [7|. The mass composition in 
these surrounding regions exhibits a strong positive corre- 
lation with environment such that more than 80% of the 
mass is in the form of resolved haloes for haloes residing in 
dense regions (for 1 + (Jt-fof ^ 5), while only ~ 20% of 
the mass is in resolved haloes for haloes residing in voids 
(for 1 -I- S7-F0F ^ 0.4). Even though only a fraction of these 
haloes enters the virial radii and contributes to the actual 
growth of the central halo, it is the environmental property 
of this mass reservoir that leads to the environmental de- 
pendence of Mnicr- For the diffuse growth component Mdif, 
however, its negative correlation with local density observed 
in Fig. |4] is not explained by the environmental dependence 
of the available diffuse mass in the reservoir outside of the 
haloes. In fact, more diffuse material is available in the reser- 
voir in denser regions (dashed curves in Fig. [7f. Why then IS 
the diffuse accretion rate Mdif lower in denser regions? We 
speculated that some of the diffuse mass may be dynami- 
cally hotter due to tidal stripping and therefore is harder 
to accrete. A careful analysis using the actual particle data 
from the simulation would be needed to investigate this issue 
further. 

The halo growth rates Mmcr and A/dif together account 
for the overall halo mass accretion history (MAH) that we 



2009 1. When combined, the opposite correlations of the two 
rates with the local density cancel each other to some extent, 
resulting in a weak environmental dependence for the for- 
mation redshift 2/ (Fig. |6|. At a fixed mass, galaxy-sized 
haloes in overdense regions on average form earlier than 
those in underdense regions, consistent with the assembly 



bias result reported in several recent papers ( Gao et al. 


2005 


Sheth & Tormen|2004 Jing et al.|2007| 


Gao & White 


2007 


Harker et al. 


2006 


Wechsler et al.|2006 


Wetzel et al. 


2007 


Wang et al.| 


2007 


Hahn et al. 20081. 


The Zf for cluster- 



sized haloes, on the other hand, show no dependence on 
environment. We note that had we ignored A^dif and only 
taken into account the mass growth due to mergers, then 
the positively-correlated S dependence of Mmcr would imply 
a later formation redshift (i.e., a smaller 2/) for haloes in 
denser regions, which is opposite to that shown in Fig. [6] 
The negatively correlated Mdif with the local density there- 
fore plays a crucial role in counterbalancing the positively 
correlated Mmcr so that the environmental dependence of 
the total rate is consistent with that of 2/. 

We have emphasized throughout the paper that the val- 
ues of Mmcr vs. Mdif depend on the mass threshold ^min used 
to defined these two components. Despite this fact, our tests 
have shown (see Fig.[8| that the enviromental trends of Mmer 
and Mdif in Figs. [4] - [t] are robust and are independent of 
the definition. 

A number of interesting questions remain to be ex- 
plored. For instance, what is the nature of the non-halo 
component Mdit? Besides tidal stripping, are there addi- 
tional physical processes controlling the different environ- 
mental dependence of Mmcr vs Mdif? How do our results for 
Mmer and Afdif extend to higher-resolution simulations, e.g., 
will the power-law slope of Mmer vs progenitor mass ratio ^ 
remain the same or steepen at smaller ^, and will this change 
be 5-dependent? How does one introduce non-Markovian 
features into the random-walk picture in the much-used EPS 
model in order to reproduce the halo growth histories as a 
function of halo mass and halo environment reported in this 
paper? We expect to shed light on some of these questions 
by using the Millennium-II simulation ( [Boylan-Kolchin et al.| 
2009 1 with ~ 1000 times better mass resolution in an up- 



commg paper. 
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